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Abstract 

T his paper considers various aspects of the initial value problem for fractional order differential 
equations The main contribution of this paper is to use the solutions to known spatially 
distributed systems to demonstrate that fractional differintegral operators require an initial condition 
term that is time- varying due to past distributed storage of information. 


Introduction 

The problem to be considered here is the initial condition problem for systems with components of 
fractional mathematical order. Specifically, the case when the fractional components are of 
mathematical order 1/2 is considered at length. This is due to the availability of a convenient 
physical system to use as reference, namely, the semi-infinite lossy line. Laplace transform 
techniques are utilized throughout the study, as they are the fundamental tool of linear system 
theory. The main contribution of this work is the demonstration that fractional differintegral 
operators require the addition of time-dependent initial condition terms in their application. The 
implication is that the solution of fractional ordinary differential equations will require that the 
problem statement include the appropriate time-dependent initial condition terms along with some 
information about the structure of these terms. It is also shown, by using a semi-infinite elastic 
bar, that the simple derivative or integral can require these terms as well. The present section 
provides some background into fractional order systems before the distributed line is addressed in 
the later sections. 


The idea of fractional integrals and derivatives has been known since the development of the 
regular calculus, with the first reference probably being associated with Leibniz m 1695 [Oldham 
and Spanier (1974), page 3]. Although not well known to most engineers, the fractional calculus 
has been considered by prominent mathematicians [Courant and Hilbert (1953)] as well as the 
"engineers" of the operational calculus [Heaviside (1893)], [Bush (1929)]. In fact many 
textbooks written before 1960 have a small section on fractional calculus [Goldman (1949)], 
[Holbrook (1966)], [Starkey (1954)], [Carslaw and Jeager (1941)], [Scott (1955)], [Mikusinski 
(1959)]. An outstanding historical survey can be found in Oldham and Spanier who also give 
what is unquestionably the most readable and complete mathematical presentation of the fractional 


l 


calculus. Other bound discussions of the area are given by Ross, McBride, and McBride and 
Roach. Unfortunately, many of the results in the fractional calculus are given in the language of 
the mathematician and are not readily accessible to the general engineering community. 

Many systems are known to display fractional order dynamics. Probably the first physical system 
to be widely recognized as demonstrating fractional behavior is the semi-infinite lossy (RC) line. 
Looking into the line, the operational impedance is 


Z(s) = 


V(s) _ 1 

I(s) Vs ' 


The current is thus equal to the half-derivative of the applied voltage (assuming time begins at t=0), 
which can be written in the time domain as, 


i(t) = 


d 1/2 v(t) 
dt 1/2 ' 


Although studied by many, Heaviside considered this system extensively by using the operational 
calculus. He also states that "there is a universe of mathematics lying in between the complete 
differentiations and integrations" [Heaviside (1893), p. 459] This specific system is studied 
further in the next section. An equivalent system is the diffusion of heat into a semi-infinite solid. 
Here the temperature at the boundary is equal to the half integral of the heat rate there. Other 
systems that are known to display fractional order dynamics are viscoelastic systems [Bagley and 
Calico (1991)], [Koeller (1984)], [Koeller (1986)], [Skaar, Michel, and Miller (1988)], [Lopez- 
Marcos (1990)]; colored noise [Mandelbrot (1967)]; electrode-electrolyte polarization [Ichise, 
Nagayanagi, and Kojima (1971)], [Sun, Onaral, and Tsao (1984)]; dielectric polarization [Sun’ 
Abdelwahab, and Onaral (1984)]; boundary layer effects in ducts [Sugimoto (1991)]; and 
electromagnetic waves [Heaviside (1893)]. As many of these systems depend upon specific 
material and chemical properties, it is expected that a wide range of fractional order behaviors are 
possible by considering different materials. 


The two most common definitions for the general fractional differintegral are probably the 

Grunwald definition and the Riemann-Liouville definition [Oldham and Spanier (1974)]. The 
Riemann-Liouville definition is given here: 


<flf(t) = 1 

[dft-a)]^ r(-q) 


f(T) 


(t - t)^ +1 


dx , q<0. 


(D 


where (a) represents the beginning of the domain of the differintegral. Here q can have non-integer 
values, and thus the name fractional differintegral. Notice that the definition is based on 
integration, and more importantly is a convergent convolution integral for q<0. For fractional 

differentiation (q>0), the usual integer n^ 1 derivative is take l of the fractional (q-n)^ 1 integral, 


d^f _ cp_ [ dq- n f 

dtq dt n IdtH ’ 


q>0 and n an integer > q . 


( 2 ) 


The basic engineering tool for analyzing linear systems i:i the Laplace transform. The Laplace 
transform of Equation ( 1 ) as given by Oldham and Spanier vith a=0 is 
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, for all q , 


(3) 



s k 


. dt^ 1 ^ 


J t=o 


where n is an integer such that n-l<q<n. The problem with this formulation is that the initial 
condition terms are fractional integrals and derivatives evaluated at time zero. Unfortunately, this 
information is not readily specified in many problems. Notice that when the initial conditions are 
assumed to be zero, this equation reduces to 


= s qi{f(t)} . ( 4 ) 

It is the nature of the initial condition terms that motivates this study. The next section addresses 
this problem through a study of the semi-infinite lossy line. 

The Semi-Infinite Lossy Line 

This section contains a classical iterated Laplace transform derivation of the operational impedance 
for a semi-infinite lossy line. Of particular importance is the inclusion of the spatial initial 
conditions in the formulation. These are necessary in order to specify the initial state completely, 
so that a proper time response can be obtained. Thus the semi-infinite lossy line becomes the 
reference case for understanding the initialization of more general fractional order differential 
equations. 

Assuming a lossy RC-line, the boundary value problem can be written in terms of the current or 
voltage variable [Karplus (1958)]. Since a semi-infinite line is considered, the measurable inputs 

or outputs are assumed to occur at the left end (x=0), while the right end (x=°°) is kept at a finite 
value. The desired result is to mathematically define the behavior of the line at x=0, that is, to 
determine both the input impedance and the initialization response. In terms of the voltage 
variables, the problem can be written as 


= a v (0,t) = Vi(t), v(oo,t) = 0, v(x,0) given, 

3x 2 


with 


i(x,t) 


1 3v(x,t) 
R dx 


(5) 


( 6 ) 


where v is the voltage, i is the current, Vi(t) is a time-dependent input variable, and Ct is a constant 
which is one over the product of resistance per unit length, R, and capacitance per unit length, C. 
A classical solution using iterated Laplace transforms is used to solve this problem. 


Taking the Laplace transform with respect to time and using s as the temporal Laplace variable 
gives 
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sV(x,s) - v(x,0) = a -- V ( x,s ) , 

dx 2 

with 


V(0,s) = V,(s), V(oo,s) = 0, I(x,s) = J-^^. ( 7) 

R dx 

Then taking the Laplace transform with respect to spatial position x, and using p as the spatial 
Laplace variable gives 

%V(p,s) - lV(p,0) = p 2 V(p,s) - pV(0,s) - . (8) 

•- dx (o,s) 

Replacing dV/dx with V*, this equation can be manipulated to give V(p,s) as 

r p 2 - ^]V(p,s) = - lV(p,0) + pV(0,s) + V*(0,s) (9) 

or 



Each term in this equation will be considered separately in order to inverse transform this with 
respect to the p variable. The first term, the transform of the initial spatial distribution, V(p,0), is 
problem-dependent. A slight rearranging of the expression gives; 



+ V-rl [p v <°.s) + v*(0,s)] . (i i) 

L r aJ 

Here the first term represents the effect of any initial spatial voltage distribution, while the second 
term represents the voltage and current present at the x=0 end of the line. The first term is now 
inverse-Laplace-transformed with respect to the p variable using convolution, and the other term is 
inverse-transformed using standard transform pairs; 
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V(x,s) = 


e + Vf (,l ' X) -i vft.O) 

S L OC 


dX - 


J o 


e-VF x 4£ v(X,0) 


dX 


Equivalently, 



V*(0,s) 




( 12 ) 


V (x,s) = 


1 


2 Vt 


- * e-V- x 
a c V a 


v(X,0) 

a 


dX 


/■* 



Ha X v(^’O) 




+ 


V(0,s) 





(13) 


Collecting like exponentials then gives 


V(x,s) 



V(0,s) + 


V*(0,s) 




e'Vf x v(X,,0) dX, 


+ 



(14) 


It should be recognized that the coefficients multiplying the two exponential functions are 
unknowns. Although the integral and either V(0,s) or V*(0,s) are given in the problem statement, 
the other condition (V* or V, respectively) at x=0 is determined as a response to these two given 
terms. Imposing the boundary condition at x = °°, allows the determination of a relationship 
between these three terms at x=0, and thus allows the impedance and initial condition response of 
the system. 


It is shown in the Appendix that in the limit as x->oo, the second term in Equation (14) goes to zero 
through application of L’ Hospital’s rule and Leibniz’s rule, and by the fact that V(°°,s)=0 from the 
problem statement. Then 


5 



= 0. (15) 


- e + y?‘ 

2 


V(0,S) + ^2^ - Lim —L= f e-Vf x v(K0) 

Va x ^~ a n Jo 


V(oo,s) 

Dropping the limit notation, by necessity we have 


<R 


V(0,s) + L_ I e-VJi V ft,0) d?, = 0 . 

VI «V« A 

Remembering that the current anywhere in the line is related to the voltage, then at x=0 

1(0,,) = . 1 dV(O-s) = . v^(0,s) 

R dx R 

and solving for voltage in terms of source current in Equation (16) gives 


(16) 


(16’) 


V(0,s) 


R 1(0, s) , i 

V£ “VI 




v(A.,0) dX 


(17) 


In evaluating the integral on the right, it is now recognized that this term is equivalent to a Laplace 
transform integral with the Laplace operator . Thus, Laplace transform tables can simplify the 
evaluation of this term as follows: 


1 


a 




e'Vf x v(5l,0) dX 



(18) 


The notation on the right side of this equation is used to indicate the evaluation procedure. First the 
initial spatial distribution v(x,0) is Laplace transformed with respect to the spatial Laplace variable 
p to give V(p,0). The integral on the left of Equation (18) is then easily evaluated by replacing the 

spatial Laplace variable p with fs/a. Using this. Equation (17) becomes 



(19) 


Notice that this contains the driving point impedance function Z(s), which is obtained by setting the 
initial condition terms to zero. 


Z(s) = 


V(0,s) 

1(0, s) 
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or, as (X = 1/RC, 


z( S) = M = ./ei 
' 1(0, s) V C Vs 


( 20 ) 


Notice from Equation (19) that the voltage V(0,s) is seen to be composed of two parts, the forced 
response due to 1(0, s) and the initial condition response due to the initial voltage distribution in the 
lossy line. 

Applying the other boundary condition of Equation (7) at x=0, allows the solution of the voltage 
transfer function V(x,s)/V(0,s) and has no effect on the terminal impedance. This is now done for 
completeness. At x=0, the boundary condition is 

V(0,s)=Vi(s). (21) 

Now, inserting Equations (21), (19), and (16’) into Equation (14), gives 


V (x,s) 



e'VoF x v(A,,0) dA, 


+ 


e'Va~ x 

2 


[V(p,0)]p=jr 
2 V t (s) -=^~ + 

1 | 

• X 

e + VoF^ v(A,,0) dA, 

I 

Qh 

a Va 1 

0 


( 22 ) 


This equation is the general expression for the voltage anywhere inside the line as a function of the 
applied voltage at the terminal Vj(s) and the initial conditions on the line. Furthermore, the current 
at any point in the line can be determined directly from Equation (22) as 


I(x,s) = - 


1 dV(x,s) 
R dx 


(7) 


We will now discuss the voltage-current behavior of the lossy line. Referring back to the 
definition of the Laplace transform from Equation (4), and applying it to Equation (19), it can be 
seen that the input-output properties (with initial conditions equal to zero) of die lossy line satisfy 


i(t) = 


1 d 1/2 v(t) 

RVa dt 1/2 


or v(t) = RVa 


d‘ 1/2 i(t) 


dt 


- 1/2 


More completely, from Equation (17) or (19), 


v(t) = RVa 


d~ 1/2 i(t) 

df 1/2 


+ <h(t) 


(23) 
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where <t>l (t) is obtained as the inverse Laplace transform of either term in Equation (18). 
Likewise, 


i(t) = 


1 d 1/2 v(t) 

RVa dt 1/2 


+ <J>2(t) 


(24) 


where <j> 2 (t) is effectively the negative integral of <t>i(t) weighted by 1/RVa. From these 
expressions it should be clear that the fractional differintegral of any function must include a time- 
dependent initial condition term. 


The Semi-Infinite Lossless Line 

This section considers the semi-infinite lossless line whose impedance is a constant or an operator 
of order zero. A classical iterated Laplace transform is again used for the derivation of the 
operational impedance for the semi-infinite lossless line. The inclusion of the spatial initial 
conditions in the formulation is of particular importance. These are necessary in order to specify 
the initial state completely, so that a proper time response can be obtained. The results will aid in 
the understanding of the initial condition terms for integer order systems in a fractional ordinary 
differential equation setting. 

Assuming a lossless LC-line, the boundary value problem can be written in terms of the current or 
voltage variable [Karplus (1958)]. Since a semi-infinite line is considered, the measurable inputs 

or outputs are assumed to occur at the left end (x=0), while the right end (x=°°) must be kept at a 
finite value. The desired result is an impedance at the left end of the line as well as the response 
due to any initial conditions on the line. In terms of the voltage variables, the problem can be 
written as 

9 2 v(x,t) _ 3 2 v(x,t) 

at 2 "lc ax 2 ’ 

v(0,t) = vi(t) , v(°°,t) = 0 

with 

di(x,t) _ _ c 3v(x,t) and 9v(x,t) _ _ L 
dx dt dx 

where v is the voltage, i is the current, Vj(t) is a time-dependent input variable, L is inductance per 
unit length, and C is capacitance per unit length. A classical solution to this problem is obtained 
using iterated Laplace transforms. 


v(x,0) and v'(x,0) given, (25) 


9i(x,t) 


(26) 
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Taking the Laplace transform with respect to time and using s as the temporal Laplace variable 
gives 


s 2 V(x,s) - sv(x,0) - ^ V ^ X — 

dt 


1 d 2 V(x,s) 
LC dx 2 


with boundary conditions 


V(0,s) = Vi(s), V(~,s) = 0, 


dV(x - s I = - L |4<M = - L [sl(x,s) - l(x,0)] . (27) 

dx v 1 dt II 

Then taking the Laplace transform with respect to spatial position, x, and using p as the spatial 
Laplace variable gives 


LCs 2 V(p,s) - LCsV(p,0) - LCV'(p,0) = p 2 V(p,s) - pV(0,s) - 


dV 


dx 


(0,s) 


(28) 


where dv(x,0)/dt is replaced with V’. Further replacing dV(0,s)/dx with V*, this equation can be 
manipulated to give V(p,s) as 

[p 2 - LCs 2 ]V(p,s) = - LCsV(p,0) - LC V'(p,0) + pV(0,s) + V*(0,s) (29) 


or 


V(p,s) 


1 

p 2 - LCs 2 


[-LCs V(p,0) -LC V'(p,0) + pV(0,s) + V*(0,s)]. 


(30) 


Each term in this equation will be considered separately in order to inverse-Laplace- transform this 
with respect to the p variable. The first term, the initial spatial distribution V(p,0), is problem- 
dependent. A slight rearranging of the expression gives 


V( P ,S) 


1 

2sVLC (p - sihc) 


1 

2sl/LC (p + sfLc) 


[-LCs V(p,0) - LC V'(p,0)] 


+ 


1 

p 2 - LCs 2 


[pV(0,s) + V*(0,s)] . 


(31) 


Here the first term represents the effect of any initial spatial distribution, while the second term 
represents the voltage and current at the x=0 end of the line. The first term is now inverse- 
Laplace-transformed with respect to the p variable using convolution, and the other terms are 
inverse-transformed using standard transform pairs; 
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V (x,s) 



e+ VLC s(x-^_ V (X,0)] dl 


2sVIC 


e +Vix s(x-X)[. V ’(^,0)] dX, 



e -VLCs(x-X)[_ v(X.,0)] dX. 


e -VEc s(x-X)[. V ’(X,0)] dX, 


+ V(0,s) coshffLC sx) + sinhfVLC sx) 

i JLL. S 


(32) 


Collecting like exponentials then gives 


V(x,s) = 


_ e 


+Vlc 


SX 


V(0,s) + m 

VLCs 


- — I e' v ^ v(A,0) dA — -1 — ( 
fLC . sfLC / 

J 0 


,-Vlc sA, V * 


v’(X.,0) dX, 


»-Vlc 


sx 


V(0,s) 


V*(0,s) 

/LCs 


+ 


fLC 


e+VEc sX v (^ 0 ) d ^ 


sfLC 


f 


e +VLC sx v '(^ 0 ) d ^ 


• (33) 


It should be recognized that the coefficients multiplying the two exponential functions are 
unknowns. Although the integrals and either V(0,s) or V*(0,s) are given in the problem statement, 
the other condition (V* or V, respectively) at x=0 is determined as a response to the given terms. 

Imposing the boundary condition at x=°°, allows the determination of a relationship between these 
three terms at x=0, and thus the impedance and initial condition response of the system. 


As was shown for Equation (14), in the limit as x-»°o, the second term in Equation (33) goes to 

zero through application of L’Hospital’s rule and Leibniz’s rule, and by the fact that V(°°,s)=0 
from the problem statement. Then 


V (°°,s) = 0 = & 


+Vlc 


s°° 


v(o,s) + XXM 
VLCs 


1 

Ilc 


6 -Vlc si V (X,,0) dX, — — 

sVLC 



c -Vlc sX 


v'(X,0) dX, 


(34) 
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and by necessity, 


V(0,s) 


, V*(0,s) 
VLCs 


1 

fLC 


0 -Vlc sX V (A,0) dA 


1 

sfLC 


e -VLC sX V '(X,,0) dA = 0 . 


Solving for the voltage in terms of the voltage gradient in Equation (35) gives 


(35) 


V(0,s) = - 


V*(0,s) , 1 

fLCs fLC 


6 -Vlc sX v (A, 0) dA 


+ 


1 

sfLC 


e-VLC sX v '(A, 0) dA 


(36) 


In evaluating the integrals, it is now recog nized that these terms are equivalent to a Laplace 

transform integral with the Laplace operator fLC S. Thus Laplace transform tables can simplify 
the evaluation of the third term as follows: 


1 

fLC 



e-VLC sX V (A,0) dA 


= T Lr[V(p,0)] p= VLCs, 

fLC 


(37) 


and similarly the fourth term. The notation on the right side of this equation is used to indicate the 
evaluation procedure. First the initial spatial distribution, v^O), is Laplace transformed with 
respect to the spatial Laplace variable p to give V(p,0). The integral on th e left of Equation (37) is 

then easily evaluated by replacing the spatial Laplace variable p with if LC S. Using this, Equation 
(36) becomes 


V(0,s) = - V *jM + — U [v (p,0)]p=VLC : 
fLC s fLC 


fLC 


[V’(p,0)]^VLC: 


(38) 


This contains the driving point transfer function (not an impedance). 


V(0,s) _ 1 

V*(0,s) fLC s 


(39) 


li 





which is obtained by setting the initial condition terms to zero. Notice from Equation (38) that the 
voltage V(0,s) is seen to be composed of two parts, the forced response due to V*(0,s) and the 
initial condition response due to the initial voltage distribution in the lossless line. Mathematically 
this is expressed in the time domain as 



where <J> 3 (t) is the time-dependent initial condition response due to the initial voltage and current 
distributions on the line, which can be obtained by inverse-Laplace-transforming the two integrals 
in Equation (36). Thus it is seen that a simple integral operator can also have time varying initial 
condition terms. 

More importantly, notice that Equation (38) contains V*(0,s). From the boundary conditions 
given in Equation (27), this term can be replaced with 


Equation (38) then becomes 



Notice from Equation (42) that the voltage V(0,s) is seen to be composed of two parts, the forced 
response due to l(0,s) and the initial condition response due to the initial voltage distribution in the 
lossless line. Considering only the first term, it can be seen that the impedance looking into this 
line is then 



which is simply a constant. Mathematically, Equation (42) s expressed in the time domain as 


dV . (X ' S) = - L [sl(0,s) - 1(0,0)] . 
dx J x =n 
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v(0,t) 



i(0,t) + <t>4(t) 


(44) 


where this (^(t) is the time-dependent initial condition response due to the initial voltage and 
current distributions on the line, which can be obtained by inverse-Laplace- transforming the three 
initial condition terms in Equation (42). Thus it is seen that even a simple constant gain operator 
can also have time-varying initial condition terms. 


The Semi-Infinite Vibrating Bar 

This section considers longitudinal vibrations in the semi-infinite elastic bar which has a transfer 
function that is a pure integral of order unity. A classical iterated Laplace transform is again used 
for the derivation of the operational impedance for the semi-infinite elastic bar. The inclusion of 
the spatial initial conditions in the formulation is of particular importance. These are necessary in 
order to specify the initial state completely, so that a proper time response can be obtained. The 
results will aid in the understanding of the initial condition terms for integer order systems in a 
fractional ordinary differential equation setting. 

For the elastic bar, the boundary value problem can be written in terms of the longitudinal 
displacement and the longitudinal force [Churchill (1958)]. Since a semi-infinite bar is considered, 
the measurable inputs or outputs are assumed to occur at the left end (x=0), while the right end 

(x=oo) must be kept stationary. The desired result is an impedance at the left end of the bar as well 
as the response due to any initial conditions on the bar. In terms of the displacement y(x,t), the 
problem can be written as 


d 2 y(x,t) _ E d 2 y(x,t) y(co,t) = 0 , y(x,0) and y’(x,0) given, (45) 

0t 2 P dx 2 


with 

_ E 3y(Q,t) _ f(t)? ( 46 ) 

8x 

where E is Young's modulus of elasticity, p is the mass per unit volume, f(t) is the force at the left 

end of the bar, and dy/dx is the strain which is proportional to the force. The force can be the 
applied input with the displacement as the output, or the displacement can be the input with the 

force as the output, as desired by the user. In what follows we will let a 2 = E/p. 

A classical solution using iterated Laplace transforms will be used to solve this problem. The 
results from this solution will provide further guidance into fractional ordinary differential 
equations. 

Taking the Laplace transform with respect to time and using s as the temporal Laplace variable 
gives 
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( 47 ) 


s 2 Y(x,s) - sy(x,0) - Mg i °) = a 2 d 2 Y(x,s) 

dt dx 2 


with 


Y(oo, s ) = 0, dY(Q ’ s) = - J- F(s) = Y*(0,s) . 
dx E 

Then taking the Laplace transform with respect to spatial position x and using p as the spatial 
Laplace variable gives 


-Ls2Y(p,s) - J-sY(p,0) - 
a 2 a 2 


-LY'(p, 0) = p 2 Y(p,s) - pY(0,s) - 
a 2 


dY 

dx. 


(0,s) 


(48) 


Replacing dY(x,s)/dx with Y* and dy(x,0)/dt by y', this equation can be manipulated to give 
Y(p,s) as 


p 2 - J- s 2 

l a 2 . 


Y(p,s) = - -L sY(p,0) - -L Y'(p,0) + pY(0,s) + Y*(0,s) 


(49) 


or 


Y(p,s) 


1 



-L s Y(p,0) - \ Y'(p,0) + pY(0,s) + Y*(0,s) . 
a 2 a 2 


(50) 


Each term in this equation will be considered separately in order to inverse-Laplace- transform this 
with respect to the p variable. The first term, the initial spatial distribution Y(p,0), is problem- 
dependent. A slight rearranging of the expression gives 


Y(p,s) 


1 


1 

IT (P ' “) 

2s | 

P + -) 

a \ r a / 

a ' 

a/J 


+ 

1 



p2-ii 
L a 2 


- -L S Y(p,o; - j- Y'(p,0) 
a 2 a 2 

[ P Y(0,s) + Y*(0,s)] . 


(51) 


Here the first term represents the effect of any initial spatial distribution, while the second term 
represents the relationship of the displacement and force at he x=0 end of the bar. The first term is 
now inverse-Laplace-transformed with respect to the p vanable using convolution, and the other 
term is inverse-transformed using standard transform pairs; 
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Y(x,s) = If e+ 


s(x-X,)/a _ y(X,,0)] dA. 


a e -s(x-?t)/a[. y(X,,0)]<a 
2 


+ I _a_ e + s ( x -*-) /a [- y'(A.,0)] dA, 
2s 

'0 


A- e - s(x-X)/a[_ y'(A,,0)] dA, 
2s 


+ Y(0,s) cosh|^j + g^’ S ^ sinh|^-j 

Collecting like exponentials then gives 


Y(x,s) = 


_ e 


+sx/a 


Y(0,s) + 


aY*(0,s) 


- a 


f 


e - sX/a y(X,0) dA, - a e* s ^ /a y'(A,,0) dX 
s 


'o 


sx/a 

+ 2 

2 


Y(0,s) - 


aY*(0,s) 


+ a 


}■* 


sWa v(A,,0) dA, + — 
s 


»+ sX/a v* 


y'(A,,0) dA, 


(52) 


(53) 


It should be recognized that the coefficients multiplying the two exponential functions are 
unknowns. Although the integrals and either Y(0,s) or Y*(0,s) are given in the problem statement, 
the other condition (Y* or Y, respectively) at x=0 is determined as a response to the given terms. 
Imposing the boundary condition at x=°o, allows the determination of a relationship between these 
three terms at x = 0, and thus allows the impedance and initial condition response of the system. 

As was shown for Equation (14), in the limit as x-»°o, the second term in Equation (53) goes to 
zero through application of L’ Hospital’s rule and Leibniz s rule, and by the fact that Y(°°,s)— 0 
from the problem statement. Then 


Y(oo, s ) = 0= e 


+ s°°/a 


Y(0,s) + 


aY*(0,s) 


- a | e" sWa y(A,,0) dX - j J tr s ^ a y'(A.,0) dA, 


(54) 
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and by necessity, 


* r f°° 

Y(0,s) 4- aY r*» s ) - a I e - si/a y(X,0) dX - 1 I e' sWa y'(A.,0) dA. = 0. (55) 

l ]„ 

Solving for the displacement in terms of the strain (or force) in Equation (55) gives 

Y(0,s) = - aY * s (0 - s) + a I e- sWa y(^,0) dX + 1 I e - *"* y'O.,0) eft . (56) 

Jo Jo 

In evaluating the integrals, it is now recognized that these terms are equivalent to a Laplace 
transform integral with the Laplace operator s/a. Thus, Laplace transform tables can simplify the 
evaluation of the third term as follows: 


a I e- ** yft.O) eft = a [YfoO)]^, (57) 

Jo 

and similarly, the fourth term. The notation on the right side of this equation is used to indicate the 
evaluation procedure. First the initial spatial distribution, y(x,0), is Laplace transformed with 
respect to the spatial Laplace variable p to give Y(p,0). The integral on the left of Equation (57) is 
then easily evaluated by replacing the spatial Laplace variable p with s/a. Using this. Equation (56) 
becomes 


Y(0,s) = ' aY * (0 ’ s) + a [YCp.O)]^ + | [YXp.O)]^ , 


(58) 


and this contains the transfer function, which is not an impedance. 


ms) _ a 
Y*(0,s) s 


(59) 


Notice from Equation (58) that the displacement Y(0,s) is seen to be composed of two parts, the 
forced response due to the applied strain and the initial condition response due to the initial 
displacement distribution in the elastic bar. Mathematically this is expressed in the time domain as 


y(0,t) = - a 


dy(0,t) 

dx 


dt + (J>5 (t) 


(60) 


where 4>s(t) is the time dependent initial condition response due to the initial displacement 
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distribution in the elastic bar, which can be obtained by inverse-Laplace-transforming the two 
integrals in Equation (56). Thus it is seen that a simple integral operator can also have time- 
varying initial condition terms. 

Likewise, if the input at the left end was displacement, the force (strain) response can also be 
determined from Equation (58); 


Y*(0,s) = - S Y(0,s) + [Y(p,0)]pw» + [Y'(p,0)Wa . 

cl 


(61) 


The strain response Y*(0,s) is then composed of two parts the forced response due to the applied 
displacement Y(0,s) and the initial condition response due to any displacement distribution on the 
bar at time zero. In the time domain, this is 


y*(0,t) = - 


1 dy(0,t) 
a dt 


+ $ 6(0 


(62) 


A gain (J) 6 (t) can be obtained from the given initial displacements. 


Summary 

This paper has demonstrated that when specific differintegral operators (q=l/2, q=l, and q=0) are 
considered as semi-infinite systems, a time-dependent term resulting from the initial spatial 
conditions should be added to the to the forced response. This is an important result with respect 
to the solution of fractional ordinary differential equations that can have rather arbitrary initial 
conditions. 

To apply this result to a specific fractional order differential equation, the initial conditions, if any, 
must be given along with some structural information on how they belong in the problem. For 
example, if the following system were given, 

d^ + d«t) + d^f(0 + f(t) = u(t)> 

dt dt 

any initial conditions must specify how they are imposed. Specifically, a physical implementation 
of this system could contain three half-integrators, or two integrators and two half integrators. The 
initial condition specification must then also specify something structural. Also, this simple 
example clearly demonstrates that further study into the meaning of minimality for fractional order 
representations is necessary. 

It should be apparent from this study that there is a clear connection between the fractional calculus 
and spatially distributed systems. Not only can spatially distributed system representations allow a 
better understanding of the fractional calculus, but the fractional calculus also allows a more 
compact representation and problem solution for some spatially distributed systems. 
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Appendix 


Equation (14) is repeated here for reference. 


V(x,s) 



V(0,s) + 



V*(0,s) 



V*(0,s) 




f 


eia l v(X,0) dX 


f 


e + Vf x v(X,0) dX 


(14) 


To get to Equation (16), which is the first term of Equation (14), set equal to zero, it is required to 

evaluate Equation (14) in the limit as x~>°°. l n this limit, the second term in Equation (14) should 
go to zero due to the exponential, however, the integral inside the bracket will diverge. That is. 



0 * °° . 


We are thus left with an indeterminate form. This can be evaluated by using L’Hospital’s rule and 
rearranging as follows 



L’Hospital’s rule says that this ratio has the same value as tie ratio of the derivatives (with respect 
to x) of the numerator and denominator. Differentiating the denominator is easy, but differentiating 
the numerator with respect to x requires the use of Leibniz s rule. Performing the differentiations 
gives 



Combining the leading constants, and applying the Leibniz’ ; rule to the numerator, gives 
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It can now be seen that the exponential terms cancel, which leaves the result 


lim 

X — >oo 


X V(x ,°) 


V(oo,0) 

2s 


The problem statement, however, requires that the boundary condition 


V(oo, t ) = 0 

be satisfied for all time. Thus it is shown that the first term of Equation (14) is equal to zero in the 
limit as x— 


19 



References 


Bagley, R.L. and Calico, R.A., "Fractional Order State Equations for the Control of 

Viscoelastically Damped Structures," J. Guid., Cont., and Dyn., vol 14, no 2, dd 304-311 
March- April 1991. 

Bush, V., Operational Circuit Analysis . Wiley, 1929. 

Carslaw, H.S. and Jeager, J.C., Operational Methods in Applied Mathematics. 2nd Oxford, 

1941 

Churchill, R.V., Operational Mathematics, McGraw-Hill, New York, 1958. 

Courant, R. and Hilbert, D., Methods of Mathematical Physics . Interscience Publishers, New 
York, 1953. ’ 

Goldman, S., Transformation Calculus and Electrical Transients. Prentice-Hall, New York, 1949. 

Heaviside, O., Electromagnetic Theory . Chelsea, New York, 1893/1971. 

Holbrook, J.G., Laplace Transforms for Electronic Engineers . 2nd, Pergamon, Oxford, 1966. 

Ichise, M., Nagayanagi, Y., and Kojima, T, "An Analog Simulation of Non-integer Order 

Transfer Functions for Analysis of Electrode Processes," J. Electroanal. Chem. vol 33 pp 
253-265, 1971. 

Karplus, W.J., Analog Simulation: Solution of Field Problems . McGraw-Hill, New York, 1958. 

Koeller, R.C., "Application of Fractional Calculus to the Theory of Viscoelasticity," J. Appl. 
Mech., vol 51, pp 299, June, 1984. 

Koeller, R.C., "Polynomial Operators, Stieltjes Convolution, and Fractional Calculus in 
Hereditary Mechanics," Acta Mechanica, vol 58, pp 251-264, 1986. 

Lopez-Marcos, J.C., "A Difference Scheme for a Nonlinear Partial Integrodifferential Equation," 
SIAM J. Num. An., vol 27, no 1, pp 20-31, Feb., 1990. 

Mandelbrot, B., "Some Noises with 1/f Spectrum, a Bridge Between Direct Current and White 
Noise," IEEE Trans. Inf. Thy., vol IT- 13, no 2, April, 1967. 

McBride, A.C., Fractional Calculus and Integral Transforms of Generalized Functions . Pitman, 
San Francisco, 1979. 

McBride, A.C. and Roach, G.F., Editors, Fractional Calci lus. Pitman, Boston, 1985. 

Mikusinski, J., Operational Calculus . Pergamon, New Yoik, 1959. 

Oldham, K.B. and Spanier, J., The Fractional Calculus . Academic Press, San Diego, 1974. 

Ross, B., Editor, Fractional Calculus and Its Applications . Springer-Verlag, Berlin, 1975. 

Scott, E.J., Transform Calculus . Harper and Bros., 1955. 


20 



Skaar, S.B., Michel, A.N., and Miller, R.K., "Stability of Viscoelastic Control Systems," IEEE 
Trans. Auto. Cont., vol 33, no 4, pp 348-357, April, 1988. 

Starkey, B.J., Laplace Transforms for Electrical Engineers , Iliffe and Sons, London, 1954. 

Sugimoto, N., "Burgers Equation with a Fractional Derivative; Hereditary Effects on Nonlinear 
Acoustic Waves," J. Fluid Mech., vol 225, pp 631-653, 1991. 

Sun, H.H., Abdelwahab, A.A., and Onaral, B., "Linear Approximation of Transfer Function with 
a Pole of Fractional Order," IEEE Trans. Auto. Cont., vol AC-29, no 5, pp 441-444, May, 
1984. 

Sun, H.H., Onaral, B., and Tsao, Y., "Application of the Positive Reality Principle to Metal 
Electrode Linear Polarization Phenomena," IEEE Trans. Biomed. Eng., vol BME-31, no 10, 
pp 664-674, Oct. 1984. 


21 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, include the time for reviewing instructions, searching existing data sources 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate for Information Operations and Reports 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0183), Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 

November 1998 


4. TITLE AND SUBTITLE 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Insights Into the Fractional Order Initial Value Problem via Semi-infinite Systems 


6. AUTHOR(S) 

Tom T. Hartley and Carl F. Lorenzo 


WU-523-22-13-00 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 11242 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM— 1998-208407 


11. SUPPLEMENTARY NOTES 


Tom T. Hartley, University of Akron, Akron, Ohio 44325-0001 and Call F. Lorenzo, NASA Lewis Research Center. 
Responsible person, Carl F. Lorenzo, organization code 5500, (216) 433-3733. 


12a. DtSTRIBUTION/A VAI LABILITY STATEMENT 

Unclassified - Unlimited 

Subject Categories: 59, 66, 31, and 70 


12b. DISTRIBUTION CODE 


Distribution: Nonstandard 


This publication is available from the NASA Center for Aerospace Information, (30 L) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

This paper considers various aspects of the initial value problem for fractional order differential equations. The main 
contribution of this paper is to use the solutions to known spatially distributed systems to demonstrate that fractional 
differintegral operators require an initial condition term that is time- varying due to past distributed storage of information. 


14. SUBJECT TERMS 



Fractional calculus; Fractional order systems; Laplace transform; Diffusion equation; 
Lossy line; Initial value problem; Semi-infinite systems 

17. SECURITY CLASSIFICATION 
OF REPORT 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

19. SECURIT f CLASSIFICATION 
OF ABSTRACT 

Unclassified 

Unclassified 

Ur classified 


15. NUMBER OF PAGES 

27 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 

















